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ABSTRACT 

We present detailed comparisons between high quahty observational colour-magnitude 
diagrams (CMDs) of open star clusters and synthetic CMDs based on MonteCarlo nu- 
merical simulations. The comparisons account for all of the main parameters which 
determine the shape of the CMD for a stellar population. For the four clusters studied, 
NGC 6819, NGC 2099 (M37), NGC 2168 (M35) and NGC 2323 (M50), we derive red- 
dening, distance, age, binary fraction, star formation rate and indicative metallicity by 
comparing the locations and density of points in the observed CMDs to the simulated 
CMDs. We estimate the uncertainties related to stellar evolution theories by adopting 
various sets of stellar models for all of the synthetic CMDs and discuss which stellar 
models provide the theoretical CMDs that best reproduce the observations. 

Key words: colour-magnitude diagrams - methods: n-body simulations - open clus- 
ters and associations: general - open cluster and associations: individual: NGC 2099, 
NGC 2168, NGC 2323, NGC 6819 



1 INTRODUCTION 



Theoretical isochrones are commonly fit to the major obser- 
vational sequences of star clusters in order to both better 
understand the underlying physics of stellar evolution and 
, to determine properties of the clusters, i.e., the age. If the 
■ metallicity, reddening and distance to the cluster are well 
constrained from independent techniques, the comparisons 
, typically involve matching the morphology of the turn-off 
and location of the red giant stars to predictions. Recently, 
the newer method of using synthetic colour-magnitude dia- 
grams to compare with observational data has proven to be 
much more informative and rewarding jAparicio et al.ll99Ci 
iTosi et al.l fl99ll : ISkillman fc GallartI '2002^ . These Monte- 
Carlo simulations allow modelling of several additional pa- 
rameters which dictate the distribution of points in the 
CMD, such as stochastic star formation (SF) processes, bi- 
nary fraction, photometric spread, main-sequence thickness, 
data incompleteness and small number statistics. Conse- 
quently, the results not only provide a measure of the prop- 
erties of the cluster, but can also constrain the star for- 
mation history (SFH) and the initial mass function (IMF). 
Furthermore, by comparing the simulations based on sev- 
eral different sets of evolutionary tracks, we can constrain 
which models use the best prescription of parameters (such 
as treatment of overshooting, mixing length, etc.). 
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Confronting the simulations with observations requires 
a large data set with accurate photometry. For this, we 
use the deep BV ph otometry pres e nted in the CFHT Open 
Star Cluster Survey llKaUrai et alJ j2001cJ) . hereafter JSKI). 
JSKI observed 19 open star clusters in our Galaxy and 
have yet published results o n the four richest clusters, NGC 
6819 jKalirai et alJ ^l2001b^ ■ hereafter JSKII), NGC 2099 
jKalirai et alJ J2001cl). hereafte r JSK III), and NGC 2168 
and NGC 2323 jKalirai et al.l ll2003h . hereafter JSKIV). 
These data were reduced and calibrated in a homogenous 
manner as described in JSKI. The resulting CMDs exhibit 
very tight main sequences showing several 'kinks' and slope 
changes which are predicted by theory. More importantly, 
the combination of very short and deep exposures, and the 
large aerial coverage of the detector (42' x 28') has allowed 
the measurement of stars from the brightest asymptotic gi- 
ant branch (AGB) and red giant branch (RGB) phases down 
to very low-mass main-sequence phases (~0.2 Mq). This al- 
lows our comparisons to yield evolutionary information over 
a wide mass range. 

The reduced data set in the CFHT Open Star Clus- 
ter Survey has been requested by, and made available to, 
several investigators for additional science rewards outside 
our goals (e.g., astrometric studies, proper motions, variable 
stars, radial velocities, brown dwarfs, blue stragglers and 
Galactic disk star distributions) . The present study comple- 
ments these efforts and analyses the four published clusters 
in a way that allows us to include them in a large homo- 
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geneous sample of open clusters aimed at studying the for- 
mation and evolution of the Galactic disk (Bragaglia 2003, 
and references therein). Galactic open clusters are indeed 
particularly well suited to this purpose, since they span a 
range of ages from a few million to several billion years and 
can be observed in various regions of the Galactic disk char- 
acterised by different star formation histories. They can be 
used to study both the present day disk structure and its 
temporal evolution (Janes & Phelps 1994, Friel 1995, Tosi 
2000, Bellazzini et al. 2003). Old open clusters offer a unique 
opportunity to trace the whole kinematical and chemical his- 
tory of our disk, if collected in populous and representative 
samples and accurately and homogeneously analysed (see 
e.g., Twarog, Ashman, & Anthony- Twarog 1997; Carraro, 
Ng, & Portinari 1998). 

Here we apply the synthetic CMD method to NGC 
6819, NGC 2099, NGC 2168 and NGC 2323 to derive their 
age, reddening, distance modulus and (approximate) metal- 
licity homogeneously to the Bragaglia (2003) cluster sample. 
The method also allows us to determine other features of 
these clusters, such as the existence (or lack thereof) of a 
significant fraction of unresolved binary systems, the origi- 
nal total mass of formed stars and the possible evaporation 
of some of the lower mass stars. 

The organisation of the paper is as follows, 321 briefiy 
summarises the data and the reduction procedures. Further 
details are given in JSKI. In f|3we present a summary of 
our main results which relate to this work from the pub- 
lished papers in the CFHT Open Star Cluster Survey. ^ 
sets up the numerical simulations and presents details on 
how the synthetic CMDs were created. Next, we compare 
the CMDs and the corresponding luminosity functions from 
the observations with the simulations on a cluster-by-cluster 
basis (JHJ- Finally, we discuss the results in JHla-nd conclude 
the study in JT] 



2 THE OBSERVATIONAL DATA 

An absolute requirement in using different data sets to make 
detailed comparisons between observations and models is 
that of homogeneity. To truly understand the variations in 
the difi'erent models parameters, we must minimize system- 
atic errors. Therefore, we desire a rich data set of open clus- 
ters, all reduced using the same approach, calibrated con- 
sistently, and with well determined incompleteness factors. 
The richest clusters from the CFIfT Open Star Cluster Sur- 
vey, NGC 6819, NGC 2099 (M37), NGC 2168 (M35) and 
NGC 2323 (M50) perfectly meet these requirements. The 
data were all imaged during an excellent three night observ- 
ing run with the Canada-France-Hawaii Telescope (CFHT) 
on 1999 October, using the CFH12K camera. The optical 
detector is a high resolution (1x10* pixels) CCD mosaic 
camera which projects to an aerial coverage of 42'x28' on 
the sky. 

Most clusters were observed with exposure times suffi- 
cient to reach the magnitude of the coolest white dwarfs and 
of the faintest M-type stars on the main sequence. We also 
complemented the deep photometry with shorter exposures 
to cover the brighter turn-off and evolved stars. The individ- 
ual images were registered and co-added (where necessary) 
using the FITS Large Images Processing Software (FLIPS) 



(J.-C. Cuillandre, 2000, Private Communication). Details of 
this procedure, as well as the pre-processing of the data are 
provided in §3 of JSKI. 

The actual photometry of the sources in the images was 
measured using Point-Spread- Function Extractor (PSFex), 
a highly automated program which will be integrated in the 
code of a future SExtractor version (Berlin & Arnouts 1995; 
E. Berlin 2000, private communication). We used a separate 
PSF on each CCD of the mosaic. Included in the PSF are 
polynomial basis functions which can map the variations of 
the PSF across the CCD. PSFex is discussed further in §4 
of JSKI. 

Data for each of the clusters were taken in the B and 
V filters, and calibrated in the Johnson photometric system 
using multiple image s of the SA-92 and SA-95 standard star 
fields jLandohll992l) . The final averaged uncertainties in the 
zero-point terms were measured to be ~0.018 magnitudes, 
uncertainties in the airmass coefficients were ~0.015 magni- 
tudes in V and ~0.007 magnitudes in B. Colour coefficients, 
as well as airmass terms, were similar to nominal CFHT ex- 
pectations. Details of the calibration of instrumental mag- 
nitudes to real magnitudes are provided in §5 and Tables 2 
and 3 of JSKI. 

The resulting, calibrated CMDs for the open clusters ex- 
hibit unprecedented long, clear main sequences ranging from 
the brightest cluster members to the limiting magnitude in 
each case (23 ^ y ^ 25). A summary of the observational 
log for each cluster (exposure time, seeing, airmass, etc..) is 
presented in Table 1. 



3 FIRST RESULTS: ISOCHRONE FITTING 

Figure 1 presents the final calibrated CMDs for each of the 
four open clusters, with the best fit theoretical isochrones 
as determined in our previous published analysis of these 
clusters (see JSKII, JSKIII, JSKIV for details). These rich 
clusters are clearly delineated from the background Galactic 
disk stars. The main sequences are very tight and show sev- 
eral features which are predicted by stellar evolutionary the- 
ory (see §5). For this, we used a new grid of model isochrones 
calculated by the group at the Rom e Observatory accord- 
ing to the prescription described in IVentura et al.l il99d) . 
These models adopt convective core-overshooting by means 
of an exponential decay of the turbulent velocity out of 
the formal convective borders fixed by Schwarzschild's cri- 
terion; a free parameter which gives the e-f olding distance 
of the exponential decay is set to ^ = 0.03 (I Ventura et alJ 
.199^) . The convective flux has been evaluated according to 
the Full Spectr um of Turbulence (FST) theory prescriptions 
iCanuto fc Mazzitelli 199^ 1. The theoretical isochrones were 
transformed into the observational plane by making use of 
the Bessel, Castelli & Pletz (1998) conversions. The lower 
main sequence (M < 0.7 Mp) was calculated by adopting 
Next Gen atmosphere models iHauschildt. AUard fc BaronI 
I1999D . For M < 0.47 M^ (or Teg < 3500 K ) the transforma- 
tions of lHauschildt. AUard fc BaronI lll999ll in B—V are not 
very reliable and so the faint end of the isochrones terminate 
at this mass. These models will be referred to as FST from 
here on. 
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Table 1. Observational Data for NGC 6819, NGC 2099 (M37), NGC 2168 (M35) and NGC 2323 (M50) 
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1. * These very short exposures were obtained at a later date. 



3.1 Published Results - NGC 6819 



The tight, very rich, main sequence and turn-off consist of 
over 2900 cluster stars to our limiting magnitude. Main- 
sequence fitting of the un-evolved cluster stars with the 
Hyades star cluster yields a distance modulus of {m—M)v 
= 12.30 ± 0.12 (d = 2500 pc), for a reddening of E{B-V) = 
0.10. These values are consistent with a theoretical isochrone 
of age 2.5 Gyrs. Detailed star counts in concentric annuli out 
to large angular radii set constraints on the cluster radius, 
R — 9'.5 rbl.'O, and clearly indicates mass segregation in the 
cluster. The global cluster mass function is found to be quite 
flat, X = —0.15 (a; — 1.35, Salpeter). A large population of 
white dwarf stars is found in the cluster and is currently 
being investigated spectroscopically. 



3.2 Published Results - NGC 2099 (M37) 

The cluster CMD shows extremely well populated and very 
tightly constrained main sequence, turn-off, and red giant 
populations. The photometry for this cluster is faint enough 
(V 24.5) to detect the end of the white dwarf cooling se- 
quence. Therefore, the white dwarfs are used as chronome- 
ters to age the cluster providing an independent age (566 
± 176 Myrs) from the main-sequence turn-off plus red gi- 
ant clump fit (520 Myrs). We also derive the reddening 
{E(B~V) = 0.21 ± 0.03) and distance ((m-M)v = 11.55 
± 0.13) to NGC 2099 by matching main-sequence features 
in the cluster to a fiducial main-sequence for the Ifyades 
ide Bruiine. Hooeerwerf fc de Zeeuwl200'lh . after correcting 
for small metallicity differences. The cluster luminosity and 
mass functions indicate some evidence for mass segregation 
within the boundary of the cluster. 
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Figure 1. Isochrone fits for each cluster are shown with corresponding ages. All models are Solar metallicity with the exception of NGC 
2168, for which Z = 0.012. For the lower main sequence, we were not able to compute non-grey atmospheres for this metallicity. We 
therefore overplot the lower main sequence from a Solar metallicity model (bottom-loft). 



3.3 Published Results - NGC 2168 (M35) 

NGC 2168 is a well studied nearby (d = 912 
i 65 Pc) young cluster (180 Myrs) showing a well 
defined main sequence. The reddening of the clus- 
ter, E{B-V) = 0.20, and the met allicity. Z = 
0.012 , were adopted from the literature iSarrazine et alJ 
I2OOOI: iBarrado v Navascues. Delivannis fc StauffeJ I2001I) . 
The cluster is found to contain ~1000 stars above our limit- 
ing magnitude with a global mass function slope very similar 
to a Salpeter value (a; — 1.29). There is mild evidence for 
mass segregation in the cluster. White dwarf stars are found 
and can potentially set very important constraints on the 
high-mass end of the white dwarf initial-final mass relation- 
ship considering they must have evolved from quite massive 
progenitors. 



3.4 Published Results - NGC 2323 (M50) 

NGC 2323 is one of the youngest clusters in the CFHT Open 
Cluster Survey (130 Myrs), and shows a main sequence ex- 
tending over 14 magnitudes in the V, B—V plane. The clus- 



ter has had very few photographic/photoelectric studies, and 
prior to our efforts, no previously known CCD analysis. This 
is surprising considering the rich stellar population (^2100 
stars above our magnitude limit) and the relatively low dis- 
tance modulus ((m— A/)o = 10.00 ± 0.17). The cluster shows 
clear evidence of mass segregation, an important result con- 
sidering the dynamical age is only 1.3 x the cluster age. 



4 SYNTHETIC COLOUR-MAGNITUDE 
DIAGRAMS 

Age, reddening and distances for the four clusters have been 
re-derived applying the synthetic CMD method (Tosi et al. 
1991) to the empirical CMDs described above. The best val- 
ues of the parameters are found by selecting the values that 
provide synthetic CMDs with morphology, number of stars 
in the various evolutionary phases and luminosity functions 
in better agreement with the empirical ones. The method 
has already been succesfuUy applied to several clusters (Bra- 
gaglia 2003 and references therein). 

The synthetic CMDs are constructed via MonteCarlo 
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extractions of (mass, age) pairs, according to an assumed 
IMF, SF law, and time interval of the SF activity. Each 
extracted synthetic star is placed in the CMD by suitable 
interpolations on the adopted stellar evolution tracks and 
adopting the Bessel et al. (1998) tables for photometric con- 
version in the Johnson-Cousins photometric system. The 
absolute magnitude is converted to a provisional apparent 
magnitude by applying (arbitrary) reddening and distance 
modulus. The synthetic stars extracted for any magnitude 
and photometric band are assigned the photometric error de- 
rived for the actual stars of the same apparent magnitude. 
Then, they are randomly retained or rejected on the basis of 
the incompleteness factors of the actual data, derived from 
extensive artificial star tests. 

Once the number of objects populating the whole syn- 
thetic CMD (or portions of it) equals that of the observed 
one, the procedure is stopped, yielding the quantitative level 
of the SF rate consistent with the observational data, for the 
prescribed IMF and SF law. To evaluate the goodness of the 
model predictions, we compare them with: the observational 
luminosity functions, the overall morphology of the CMD, 
the stellar magnitude and colour distributions, the number 
of objects in particular phases (e.g., on the red giant branch, 
in the clump, on the blue loops, etc.). A model can be con- 
sidered satisfactory only if it reproduces all of the features of 
the empirical CMDs and luminosity functions. Given the un- 
certainties affecting both the photometry and the theoretical 
parameters (stellar evolution tracks included), the method 
cannot provide strictly unique results; however, it allows us 
to significantly reduce the range of acceptable parameters. 

In this way we derive the age, reddening, distance of the 
cluster, and indicate the metallicity of the stellar evolution 
models in better agreement with the data. In principle, the 
latter should be indicative of the cluster metallicity, but this 
depends significantly on some of the stellar model assump- 
tions, such as opacities. 

To estimate if, and how many, unresolved binary sys- 
tems could be present in the cluster, the synthetic CMDs 
have been computed both assuming that all the cluster stars 
are single objects and that an (arbitrary) fraction of them 
are members of binary systems with random mass ratio. 
Again, the comparison of the resulting main-sequence mor- 
phology and thickness with the observed ones allows us to 
derive information on the most likely fraction of binaries. 

We also infer the astration mass of the cluster (i.e., the 
total mass that went in all of the stars formed in the cluster), 
according to the adopted IMF. The method here also allows 
us to derive the best IMF as long as the observed mass range 
of main-sequence stars is sufficiently large and the main se- 
quence is sufficiently tight. However, the high background 
contamination affecting the systems examined here and the 
lack of sufficiently reliable external blank fields, make the 
derivation too unsafe to discuss it here. Other factors (see 
below), such as the actual size of the photometric errors 
also prevents us from comparing the IMF in any significant 
detail. 



5 OBSERVATIONS VS THEORY - 

COLOUR-MAGNITUDE DIAGRAMS AND 
LUMINOSITY FUNCTIONS 

To test the effect of different input physics on the derived 
age, reddening and distance modulus, we have run the simu- 
lations with three different sets of stellar evolutionary tracks. 
The adopted sets were chosen because they assume differ- 
ent prescriptions for the treatment of convection and range 
from no overshooting to rather high overshooting from con- 
vective regions. Despite these differences, they are all able 
to well reproduce the observed CMDs of both star clus- 
ters and nearby galaxies. They are thus suitable to eval- 
uate the intrinsic uncertainties still re lated to stellar evolu- 
tion models. We use the FST tracks of lVentura et al.l lll998h 
with high ((^ = 0.03) and moderate (^ = 0.02) overshooting ; 
the BBC94 tracks o f the Padova group iBressan et al.llT993l : 
iFagotto et al."l994'), wit h over shooting; and the FRANEC 
tracks (Domingucz ot al] Il999|) with no overshooting. For 
those cases where the metallicities are not well constrained, 
we allow different values in our comparisons. Incompleteness 
factors and photometric errors are taken from the original 
papers in the CFHT Open Star Cluster Survey series (JSKII, 
JSKIII and JSKIV) and folded into the numerical simula- 
tions. 

All models assume that the star formation activity has 
lasted 5 Myr (i.e., approximately an instantaneous burst) 
and that the stars were formed following a single slope {x — 
1.35) Salpeter's IMF over the whole mass range covered by 
the adopted tracks (0.6 - 100 Mq). 

Figures 2, 4, 6 and 8 show for each cluster the synthetic 
models in better agreement with the data for each adopted 
set of stellar tracks, as well as the cluster observational CMD 
on the same scale. In all of these figures the CMDs in the 
top row are the synthetic ones, those in the second row are 
the same synthetic ones but overlapped with the stars of 
the equal area fiducial blank field, and the bottom CMD is 
the empirical one. Figures 3, 5, 7 and 9 compare the empir- 
ical luminosity functions with those predicted by the best 
synthetic models. Ultimately, the best models have been se- 
lected by visually comparing features between the synthetic 
models and the observations, such as the location and den- 
sity of points on the CMD. However, we have also attempted 
to confirm our results by using both the two-dimensional 
Kolmogorov-Smirnov test on the CMD and tests on the 
luminosity and colour functions. In the cases of well con- 
strained field contamination, such as behind NGC 6819 and 
NGC 2099, the tests confirm our choices of best synthetic 
models. For the other two clusters, the tests do not help 
discriminate between best choice of models and we strictly 
use conclusions based on the visual location of points. Also, 
these statistical tests are particularly not well suited for 
these comparisons because of the insignificant weight given 
to the small number of stars in key evolutionary stages which 
must be fit by the synthetic models (such as the red giant 
clump). 



5.1 NGC 6819 

NGC 6819 is the richest cluster in this work. The cluster field 
contains 10510 bona fide stars (i.e., with stellarity parameter 
^ 0.5) while the blank field of same area contains 8504 bona 
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Figure 2. Synthetic CMDs with Z = 0.02 and a 20% binary fraction (random mass ratio) are shown for NGC 6819. The BBC94 tracks 
are found to provide the best agreement with the observations (bottom panel). The data has been truncated at the faintest magnitudes 
(see Figure 1 for full CMD). See HS.ll for more information. 



fide stars. The synthetic CMDs are therefore computed with 
2000 objects (a fraction of which are assumed to be members 
of binary systems). 

All of the tracks used to simulate the NGC 6819 CMD 
are Solar metallicity models as the cluster metallicity has 
been const rained to this value th rough the high resolution 
analysis of lBraeaglia et all ||200 j) . Figure 2 shows that the 
three sets of models provide slightly different ages (and 
hence, reddening and distances), which strictly depend on 
the assumed amount of overshooting: the higher the over- 
shooting, the older the age. This is of course due to higher 
overshooting models having higher stellar luminosities, thus 
allowing lower mass stars to reach the desired luminosity. 
The variations in the ages seen in the models also stress the 
importance of deriving uncertainties on cluster parameters 
by using different sets of tracks. 

The FST models with C = 0.03 (left) lead to param- 
eters virtually identical to those found in JSKII. A better 
estimate of the age is found to be 2.30 ± 0.15 Gyrs, distance 
modulus (m-M)o = 12.0, and reddening E{B-V) = 0.10. 
The FRANEC models (right) suggest younger ages, due to 
their lack of overshooting, 1.6-1.8 Gyrs, and therefore need 



a slightly higher reddening {E{B—V) — 0.15) to reproduce 
the colours. The BBC94 models (middle), with intermedi- 
ate overshooting, favor an intermediate age of 2.0 Gyrs, with 
E{B-V) = 0.12 and (m-M)o = 12.2. 

None of these models are found to perfectly reproduce 
all of the observed features in the NGC 6819 CMD. For 
example, the FST models predict a main-sequence slope 
steeper than the observed one and the two diverge at the 
faint end. These models also tend to place the red giant 
clump excessively towards the red, probably due to the fact 
that the models don't include mass loss. The FRANEC 
models better reproduce the main-sequence distribution but 
present a turn-off with a hook that is too pronounced, a 
likely consequence of the young age. Given the high field con- 
tamination, however, this feature is not particularly evident 
in the combined CMDs of the second row. A smaller clus- 
ter core CMD shows the hook prominently. The FRANEC 
models also overpopulate the lower red giant branch. This 
is visible even despite the fact that this region of the CMD 
is plenty of contaminating objects. The BBC94 models pro- 
vide the best fit to the data CMD, despite a main-sequence 
slope which is slightly too steep. Most of these discrepancies 
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between the observed features and the synthetic CMDs re- 
suh from problems in the stellar models. As pointed out by 
Andreuzzi et al. (2004) studying the open cluster NGC 6939 
with the same method, assuming a reddening dependence on 
the stellar colour as adopted by Twarog et al. (1997) makes 
the synthetic main sequences even steeper and, hence, more 
inconsistent with the data. 

Another evident result that we find for all of the clus- 
ters studied here, is that the observed main-sequence spread 
cannot be reproduced if one assumes photometric errors as 
small as those listed in our photometry catalogue. To test 
the actual size of the photometric errors we revert back to 
the artificial star tests presented in JSKII. A plot of the Vin - 
Vout distribution tells us that, in fact, the errors quoted agree 
with the artificial star tests. The synthetic CMDs however, 
require that either the errors are larger, at least 0.01 mag- 
nitude even in the brighter bins, or there is some physical 
cause for the observed spread (i.e., differential reddening or 
metallicity spread). A spread due to binaries is modelled 
and therefore ruled out as a cause for this discrepency. The 
apparent spread could be caused by small offsets in the cal- 
ibration of the different CCDs of the mosaic, but internal 
metallicity spread in the cluster or small differential redden- 
ing cannot be excluded. We are currently investigating this 
possibility in another project (R. S. French et al. 2003, pri- 
vate communication) . For now, we have increased the errors 
to 0.01 magnitude for stars brighter than V = 19 and dou- 
bled them for fainter bins (cr = 0.01 at V = 21, a = 0.04 at V 
= 23). Empirically, this is the minimum size required to ob- 
tain a main-sequence spread comparable with the observed 
one. 

All of the shown cases in Figure 2 assume that 20% of 
the stars are in binary systems with random mass ratio. We 
also tested simulations with 30% binaries and found minor 
differences. Given the high field contamination, it is difficult 
to understand which fraction and assumption of the mass ra- 
tio is the best case. It could entirely be true that the fraction 
is actually smaller if the mass ratio is more skewed towards 
equal masses. Models without binaries are not shown as they 
amplify the problem of the main-sequence width (the pre- 
dicted one is too tight) and do not correctly reproduce the 
stellar distribution on the red side of the main sequence. 

Figure 3 shows the empirical luminosity function of the 
10510 cluster field stars (empty circles) compared to the lu- 
minosity functions resulting from the sum of the 2000 syn- 
thetic stars of Figure 2 with the 8504 stars of the blank field. 



The agreement is excellent, except for the slight underesti- 
mate of the clump stars of the FST models. 

A clear result of our study is that all of the synthetic 
CMDs for NGC 6819 show overpopulated faint main se- 
quences {V > 20), both for single and binary stars. This is 
a clear signature of low mass stellar evaporation (or at least 
segregation outwards of our cluster boundary, 9^5). As we 
showed in JSKII, the mass function of this cluster steadily 
marches to flatter values, and becomes inverted in the out- 
ermost annuli. Therefore, the observations provide the an- 
swer to this overabundance. However the effect is not seen in 
the luminosity function comparisons, which are in excellent 
agreement with the observed one. This suggests either that 
the faint magnitude luminosity function is dominated by 
field contamination or that the apparent over-concentration 
of the synthetic main sequences is actually an artifact of the 
residual underestimate of the photometric errors. 

Assuming a single slope Salpeter IMF, the mass of gas 
that formed the stars of the surveyed area of NGC 6819 be- 
tween 0.6 and 100 Mq is ~ 4.2 x IO^Mq with the BBC94 
models, 4.0 x 10^ Mq and 4.9 x 10^ M© with the FST and the 
FRANEC models, respectively. Hence, different stellar mod- 
els do not provide the same answer for a global parameter 
such as the star formation either. The differences are again 
related to the adopted input physics and to the amount of 
assumed overshooting in the stellar models. Less overshoot- 
ing implies higher masses for the observed living stars, i.e. 
shorter lifetimes, which, convolved with the negative slope 
of the IMF, leads to a higher total astrated mass. 



5.2 NGC 2099 

The cluster field contains 12194 bona fide stars, while the 
blank field extrapolated to the cluster area contains 10576 
stars. Hence, the synthetic CMDs have been computed with 
1618 objects. 

In JSKIII, we assumed the cluster metallicity of NGC 
2099 to be Solar , based on the isochrone fitting results of 
iMermilliod et all lll996l) as well as our own best gue ss from 
the fit. A more recent studv. lNilakshi fc Sagail ||200 jl . prefer 
a subsolar Z = 0.008 metallicity, again based on isochrone 
fitting. Therefore, one goal is to constrain the metallicity of 
the cluster based on synthetic CMD co mparisons. So, f or this 
cluster we use the the FST tracks of I Ventura et alJ lll998h 
with C = 0.03, Z = 0.02 and Z = 0.006, the BBC94 tracks of 
the Padova group with Z = 0. 02 (Brcssan ct al. 1993) and 
Z = 0.008 llFagotto et al.lll994^ . and the FRANEC tracks 
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Figure 4. Synthetic CMDs with subsolar metallicities and a 20% binary fraction (random mass ratio), are shown for NGC 2099. The 
BBC94 tracks are found to provide the best agreement with the observations (bottom panel). The data has been truncated at the faintest 
magnitudes (see Figure 1 for full CMD). See i|5.2l for more information. 



JPominguez et alJll999l) with Z = 0.02, Z = 0.01 and Z = 
0.006. 

For the FST, the synthetic sequences using Solar metal- 
licity models look impressively consistent with the observed 
CMD in all the minimal details. However, the turn-off and 
clump shape are not very good: the sequence after the turn- 
off is too long and curved and the clump is a bit short and 
faint. It is also very difficult to single out the best age, as 
anything between 520 and 700 Myrs doesn't significantly 
change the important CMD features (clump morphology and 
luminosity, shape of turn-off, etc.), all looking reasonable. 
The clump becomes too faint below 500 Myrs, indicating 
that this is a lower age limit and too blue beyond 800 Myrs 
indicating an upper limit. The best fit resulting parame- 
ters for these tracks are therefore age — 620 ± 60 Myrs, 
E{B-V) = 0.21 ± 0.02, and (m-M)o = 10.70 ± 0.20. The 
BBC94 Z — 0.02 tracks excellently reproduce all of the main 
sequence and clump features, but have long post turn-off se- 
quences, as we saw in the FST models. For Solar metallicity, 
these tracks are in better agreement with the data than the 
FST models. The age is also better constrained, thanks to 
a higher sensitivity of the various features: age — 590 ± 10 



Myrs, E{B-V) = 0.18 ± 0.01, and (m-M)o = 10.50 ± 
0.10. The FRANEC Z = 0.02 tracks provide the best agree- 
ment with the observed upper main sequence and turn-off 
features. However, the resulting clump is too vertical and 
overpopulated. Part of the latter problem is due to the fact 
that their main sequence is not as deep as the others, be- 
cause the minimum mass available to us is 0.7 M©, instead 
of 0.6 Mq as for the FST and Padova models. The code 
therefore can not distribute the required 1618 faint stars in 
these bins and, inevitably, puts a larger number of brighter 
stars according to a Salpeter IMF. The resulting best fit pa- 
rameters in the FRANEC models are age — 590 ± 10 Myrs, 
E{B-V) = 0.16 ± 0.02, and (m-A/)o = 10.31 ± 0.10. 

Figure 4 presents the lower metallicity synthetic CMDs 
for NGC 2099. For the Z = 0.006 FST models (left), the 
length of the sequence just after the turn-off shortens and 
becomes more vertical, thus better reproducing the observed 
morphology. The clump morphology also improves, thanks 
to a slightly larger extension. However, the agreement in the 
main-sequence morphology is no longer as excellent as the 
higher metallicity case, but still quite good. It may well be 
that intermediate metallicity tracks would work even better. 
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Figure 5. The empirical NGC 2099 luminosity function (open circles) is compared with the predicted luminosity function from the 
numerical simulations (solid line). See i|5.2l for a discussion of these results. 



Therefore, overall the synthetic CMDs from the FST tracks 
prefer the lower metallicity, with parameters: age = 520 ± 
40 Myrs, E{B-V) = 0.36 ± 0.01, {m-M)o = 10.40 ± 0.10. 
The BBC94 lower metallicity tracks (middle), Z = 0.008, 
are the best models for this cluster. They reproduce all of 
the main-sequence features and also have a better shaped 
turn-off than with the Solar metallicity models. A minor 
inconsistency is seen in the location of the clump, which is 
slightly too vertical, but this is a minor detail. The age is 
less well defined, with final best parameters age = 430 ± 
30 Myrs, E{B-V) = 0.36 ± 0.03, (m-M)o = 10.50. The 
lower metallicity FRANEC tracks (right) can be ruled out 
as both the Z — 0.01 and Z — 0.006 models show extremely 
extended clumps (real blue loops) which are inconsistent 
with the observed clump for any reasonable age. Also, the 
lifetimes in that phase favor the stars to reside on the blue 
side of the clump (vertical), adding further inconsistency to 
the colour and shape. 

We have also perfomed several tests on the fraction of 
binary stars in NGC 2099. At first glance, the empirical 
CMD doesn't show any evidence of binaries, but we have 
found a posteriori that, without including them in the syn- 
thetic CMDs, both the turn-off morphology and the colour 
distribution on the red side of the main sequence are not 
well reproduced. Our best guess is that binaries are around 
15-20% of the cluster stars and all of the results that we have 
discussed above assume this fraction (20%) with a random 
mass ratio. 

Figure 5 shows the empirical luminosity function of the 
12194 cluster field stars (empty circles) compared to the lu- 
minosity functions resulting from the sum of the 1618 syn- 
thetic stars of Figure 4 (subsolar metallicity) with the 10576 
stars of the extrapolated blank field. The predicted lumi- 
nosity functions are consistent with the data, but the agree- 
ment at the brightest magnitudes is not exceptional. The 
very good match at the faintest magnitudes shows that the 
extrapolation of the number of stars observed in the blank 
field to the larger area covered by the cluster is correct. 

The general fits of the data to synthetic CMDs favor 
the lower metallicity models for NGC 2099, thus requiring a 
higher reddening. These are found to better reproduce the 
shape of the upper main sequence and turn-off. Notice that 
age and modulus derived here do not coincide with those 
derived JSKIII with the same stellar tracks. The difference 
in age is probably due to the fact that with the isochrone 
fitting method one simply fits a curve to the points, whereas 
with the synthetic CMDs one has to reproduce the points 



themselves, including morphology and quantity. The differ- 
ence in modulus is, presumably, simply the consequence of 
the age difference. 

Assuming a single slope Salpeter IMF, the mass of gas 
that formed the stars of the surveyed area of NGC 2099 
between 0.6 and 100 M@ is ~ 2.6 x 10^ Mq with the BBC94 
models, 2.5 x 10^ Mq and 3.1 x 10^ with the FST and 
the FRANEC models, respectively. 



5.3 NGC 2168 

The cluster field contains 9298 bona fide stars. Our previ- 
ous work on NGC 2168 (JSKIV) suffered due to the large 
size of the cluster (R > 20'). A very small blank field was 
constructed from the outer edges of the CFH12K CCD mo- 
saic. In fact, some cluster members surely resided within 
the blank field, and therefore the cluster luminosity function 
could have been underestimated. Using synthetic CMDs, we 
find that we can only reproduce the observed stellar distri- 
bution in the CMD and the luminosity function if we use 
500-600 stars down to V = 19. Using the blank field from 
the outer edges of the CCD subtracts off too many stars and 
forces a poor agreement between the theory and the obser- 
vations. This confirms our initial suspicion that the cluster is 
larger than our aerial coverage. The synthetic CMDs shown 
here contain 550 stars brighter than V = 19. 

For the synthetic CMD comparisons, we again use both 
Solar and subsolar metallicity models from each of FST, 
BBC94, and FRANEC. We immediately find that the sub- 
solar metallicity models provide a much nicer agreement to 
the observations than Solar metallicity models. The Solar 
metallicity models not only predict an insufficient number 
of main-sequence stars just below the turn-off, they also con- 
tain very vertical upper main sequences. Besides, for ages as 
old as ~180 Myrs the turn-offs are slightly too faint, while 
for younger ages, <120 Myrs, the turn-off has the proper 
luminosity but is severely underpopulated when compared 
to the observations. Flattening the IMF did not help the 
overall fit. 

Figure 6 shows the best synthetic CMDs for NGC 2168. 
For the FST tracks (left), the models with maximum over- 
shooting ((" — 0.03) and low metallicity (Z = 0.006) repro- 
duce the data very nicely. These require a best fit distance 
modulus slightly smaller {{m—M)o = 9.6) and a reddening 
slightly larger {E{B-V) = 0.22) than found in JSKIV. The 
differences can be most likely attributed to the lower metal- 
licity (JSKIV uses Z = 0.012) that we chose in these com- 
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Figure 6. Synthetic CMDs with subsolar metallicities and a 30% binary fraction (random mass ratio), are shown for NGC 2168. The 
BBC94 tracks are found to provide the best agreement with the observations (bottom panel). The data has been truncated at the faintest 
magnitudes (sec Figure 1 for full CMD). See i|5.3l for more information. 



parisons. The best age for this model is 200 Myrs, slightly 
older than what JSKIV found (180 Myrs). The absence of 
a population of post main-sequence stars coupled with the 
insensitivity of the turn-off luminosity and morphology with 
age in these models allows for a nice agreement between the 
theory and data with ages down to 120 Myrs. 

The agreement is found to be not as good if we use the 
same metallicity FST models with slightly lower overshoot- 
ing, ^ = 0.02. At an age of 180 Myrs the turn-off is too faint, 
however, with younger ages (150 Myrs) we overpopulate the 
post main-sequence stars (five are seen). At 120 Myrs, the 
turn-off is too bright, hence this age is too young. 

The BBC94 models for NGC 2168 are shown in the 
middle panel of Figure 6. The Z = 0.08 tracks are found to 
excellently reproduce both the CMD morphology and the 
number of stars in the different phases. The main-sequence 
curves are best fitted using E(B—V) = 0.20 and {m—M)o = 
9.8. The best age estimate is found to be 180 Myrs. An age 
of 150 Myrs doesn't provide any post main-sequence stars 
whereas an age of 200 Myrs is too old as it predicts too many 
clump stars. The FRANEC models (right) predict too many 
post main-sequence stars for both adopted metallicities, Z = 



0.01 and Z = 0.02. As summarised in JSKIV, many studies 
of NGC 2168 exist in the literature with none finding more 
than a few clump stars. The FRANEC models would require 
a very young age (110 Myrs) to fit the post main-sequence 
phases and this would in turn force a poor agreement in the 
turn-off region. The best fit model is, arguably, the one with 
age = 110 Myrs, E{B-V) = 0.22, and (m-M)o = 9.8. 

Our best estimate for the NGC 2168 binary fraction is 
found to be ~30%. This fraction puts a consistent number 
of stars, and sparse regions, on the red side of the cluster 
main-sequence in the synthetic CMDs. All CMDs shown in 
Figure 6 contain this binary fraction. 

Comparing the resulting best synthetic CMDs with the 
observed one, it is clear that the predicted main sequence is 
too tight and slightly underpopulated at faint magnitudes 
(V > 15). As addressed earlier in t|5.1l the tightness is a 
question of photometric errors, possible differential redden- 
ing, or metallicity spreads within the cluster. The underpop- 
ulation in this case is most likely related to the uncertain 
method we have used to infer the number of cluster mem- 
bers. Many tests that we have done indicate that we can- 
not attribute more stars to the cluster without getting too 
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Figure 7. The empirical NGC 2168 luminosity function (open circles) is compared with the predicted luminosity function from the 
numerical simulations (solid line). See i|5.3l for a discussion of these results. 



many bright stars, both on the main sequence and in subse- 
quent evolutionary phases. We also note here that V — 19 
is the limiting magnitude of the numerical simulations for 
this cluster. This i s related to the ava ilability of photometric 
conversion tables iBessell et al .119981) for the adopted metal- 
licities. With that said, it is remarkable how well the overall 
shape of the upper main sequence is reproduced. The 'kinks' 
and slope changes all fall in the correct place as dictated by 
the observations. This allows us to accurately determine the 
best fit parameters for NGC 2168. Interestingly, the best fit 
synthetic CMD for this cluster, the BBC94 model, provides 
an identical reddening, distance, and age to what we found 
using the FST isochrone in JSKIV. The FST synthetic CMD 
derived parameters are also in excellent agreement with our 
observations considering that the metallicities are slightly 
different. 

Figure 7 shows the empirical luminosity function of the 
9298 cluster field stars (empty circles) compared to the lu- 
minosity functions resulting from the sum of the synthetic 
stars of Figure 6 with the 8748 stars resulting by randomly 
multiplying the 1150 objects of the supposed blank field for 
an area normalisation factor. Clearly, a major side effect of 
the lack of an appropriate external field is that the com- 
parison of the luminosity functions becomes only a matter 
of self-consistency and cannot be used for any consideration 
on mass segregation/evaporation. As mentioned above, the 
luminosity functions agree with the data only with 500-600 
cluster members. One may think to improve the situation 
by steepening significantly the IMF (thus allowing for more 
stars in the fainter main sequence, without adding bright 
ones), but we consider it highly dangerous to modeify the 
IMF without a good field subtraction. We stress that obser- 
vations of a deep BV field offset from the centre of NGC 
2168 by at least one degree would be very useful. 

Assuming a single slope Salpeter IMF, the mass of gas 
that formed the stars of the surveyed area of NGC 2168 
between 0.6 and 100 Mq is ~ 8.5 x 10^ Mq with the BBC94 
models, 7.8 x IO^Mq and 9.6 x 10^ Af© with the FST and 
the FRANEC models, respectively. 



5.4 NGC 2323 

The cluster field contains 11115 bona fide stars. Taking into 
account that 1) the blank field contains 5142 stars but covers 
an area 1.8 times smaller than the cluster field, 2) the mass 
of the available stellar models don't allow synthetic stars 
fainter that V ~ 18, the synthetic CMDs have been initially 



computed with 565 objects brighter than V = 18. A poste- 
riori, we have however found that in this way we predict lu- 
minosity functions clearly overestimating the bright portion 
(where the cluster dominates) and underestimating the faint 
portion (where field contamination dominates). We empiri- 
cally find that good luminosity functions are obtained only 
if the cluster contains 285 members brighter than V = 18, 
suggesting that the previously used blank field underesti- 
mated the actual contamination. Again observations of an 
external field of appropriate size would be crucial for a safer 
derivation of some of the cluster properties. 

For the youngest cluster in this work, NGC 2323, the 
metallicity is known to be Solar and in fact, we find that the 
Solar metallicity tracks provide the best agreement to the 
data in this work. The most difficult parameter to fit in NGC 
2323 is the cluster age. Although ages can be constrained 
from the turn-off morphology, the luminosity of the red gi- 
ant clump is an alternative tool. NGC 2323 only shows two 
bright, red objects which are potential post main-sequence 
candidates. However, since these two objects do not have 
any membership probability from the literature, we can not 
trust them as cluster members. Furthermore, the turn-off 
itself is not heavily populated preventing an accurate deter- 
mination of its luminosity. As we'll see below, we find that 
models with a range of young ages are equally consistent 
with the observed distribution. 

Figure 8 presents our best synthetic CMDs for NGC 
2323. All of the shown cases are for Solar metallicity and 
assume a 30% binary fraction with random mass ratio. 

The FST models (left) are found to perfectly reproduce 
the majority of the observed features. The agreement is also 
found to be independent of the adopted overshooting — 
0.02 or C = 0.03). Affhough the age is slightly higher (by ~20 
Myrs) in the higher overshooting case, the reddening, dis- 
tance modulus and goodness of fit are the same. Our best pa- 
rameters are therefore age = 130 Myrs, E{B—V) = 0.22, and 
(m-M)o = 10.0 for C = 0.03. However, ages between 120 
and 180 Myrs are equally likely. We have chosen 130 Myrs as 
this age predicts a post main-sequence star where a real star 
(member or not) is indeed located on the CMD. Therefore, 
the results from this work for the FST models are identical 
to the parameters found in JSKIV. The BBC94 stellar tracks 
also perfectly reproduce features on the CMD of this cluster. 
For E{B—V) — 0.22, the best fit distance modulus is found 
to be {m—M)o = 10.2. In practice, all ages between 80 and 
150 Myrs fit the data equally well. We've again chosen 120 
Myrs as the best age as it predicts a post main-sequence star 
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Figure 8. Synthetic CMDs with Z = 0.02 and a 30% binary fraction (random mass ratio), are shown for NGC 2323. The BBC94 tracks 
are found to provide the best agreement with the observations (bottom panel). The data has been truncated at the faintest magnitudes 
(see Figure 1 for full CMD). See H5.4l for more information. 



in the correct location as dictated by the potential red giant 
in the data. The FRANEC models (right) predict a larger 
number of post main-sequence stars than the FST and BBC 
models. This is therefore inconsistent with the observational 
data. For Solar metallicity, the best compromise between the 
number of post main-sequence stars, the turn-off luminos- 
ity, and the shape of the upper main sequence is found for 
an age of 110 Myrs. The turn-off is still too faint, however, 
younger ages produce very vertical and blue upper main se- 
quences. The best fit reddening is E{B—V) — 0.22 and the 
distance is {m—M)o = 10.0. Since the results from these 
latter comparisons were not very satisfactory, we also tried 
a lower metallicity, Z = 0.01. The results are very similar, 
with an overproduction of post main-sequence stars and a 
worse shape to the upper main sequence. 

Figure 9 shows the empirical luminosity function of the 
11115 cluster field stars (empty circles) compared to the lu- 
minosity functions resulting from the sum of the synthetic 
stars of Figure 8 with the 10576 stars of the extrapolated 
blank field. As already discussed for NGC 2168, the lack 
of an appropriate estimate of the back/foreground stars 
doesn't allow for a critical analysis of the luminosity func- 



tion. We can see, however, from Figure 8, that within the 
uncertainties, the selected models reproduce quite well the 
observed luminosity function. 

Assuming a single slope Salpeter IMF, the mass of gas 
that formed the stars of the surveyed area of NGC 2323 
between 0.6 and 100 Mq is ~ 4.8 x 10^ Mq with the BBC94 
models, 4.3 x 10^ M© and 5.0 x 10^ with the FST and 
the FRANEC models, respectively. 



6 DISCUSSION 

Table 2 summarises the results described in the previous 
section and includes other information about the four clus- 
ters studied here. We recall that stellar models that include 
overshooting from convective regions imply older ages than 
models without overshooting. This explains why our ages 
resulting from the FST tracks with ^=0.03 are systemati- 
cally older than those obtained with the BBC94 models and 
much older than those obtained with the FRANEC ones. 
It is interesting to notice that while the FST models with 
higher overshooting reproduce the data better than those 
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Figure 9. The empirical NGC 2323 luminosity function (open circles) is compared with the predicted luminosity function from the 
numerical simulations (solid line). See i|5.4l for a discussion of these results. 



with (^—0.02, the BBC94 models appear to be those in better 
overaU agreement, in spite of having a formal overshooting 
somewhat lower than the FST (^=0.03 models. This indi- 
cates that what counts in the predictions of stellar evolu- 
tion models is the combination of the many parameters and 
input physics. We find that in general the BBC94 models 
better reproduce both the morphological features and the 
numbers of stars in the various evolutionary phases of the 
observed clusters. 

The excellent delineation of the various main-sequence 
curves and kinks allows us to put fairly stringent limits on 
both the reddening and the distance modulus of each clus- 
ter, independent of the age. Table 2 shows in fact that the 
variation of the values of these two quantities resulting from 
different sets of stellar tracks is quite small. This clearly 
helps to better define the age as well. 

As discussed for the case of NGC 6819, the synthetic 
main sequences of all the four clusters are tighter than the 
observed main sequences when the formal photometric er- 
rors of the survey catalogues are assumed. To reproduce the 
observed width, the synthetic CMDs require either photo- 
metric errors of at least 0.01 magnitude even in the brighter 
bins or some other factors such as an internal metallicity 
spread or differential reddening. 

The synthetic CMD tests clearly indicate that NGC 
2168 is larger than previously thought and that our outer 
CCD blank field is not adequate. A similar problem, al- 
though to a lesser degree, may also exist in the NGC 2323 
data set. Future deep photometry of these clusters should 
include separate blank field observations. Indeed, in spite 
of the great depth and excellent accuracy of our photome- 
try, which allows us to measure the faintest, least massive 
stars of each cluster, the lack of appropriate decontamina- 
tion fields has prevented us from safely deriving their present 
mass functions and, consequently, their IMF, mass segrega- 
tion, and evaporation status. 

The synthetic CMDs have allowed us to infer the as- 
trated mass of each cluster. We find that although this does 
depend on the adopted set of stellar evolutionary tracks, it 
only varies within a factor 1.2 from the minimum to the 
maximum value obtained with the stellar models considered 
here. We have chosen to simply apply Salpeter's single slope 
IMF for this exercise. We note, however, that increasing evi- 
dence suggests that the slope is most likely different (steeper 
for more massive stars, flatter for less massive ones, e.g., 
Kroupa 2002) beyond the range of masses, 2 — IOMq, orig- 
inally analysed by Salpeter (1955). Although this simplistic 



assumption introduces a systematic error in our estimates, 
it is not a large concern as we are not providing absolute 
values of star formation rates but rather just comparing rel- 
ative rates for the four clusters. 

Given the angular size of the observed field of view for 
each cluster (see Table 2), the actual area surveyed obviously 
depends on the adopted distance to the cluster. We have as- 
sumed both the instrinsic distance modulus and the astrated 
mass resulting from the best synthetic CMD comparisons, 
the BBC94 models. This corresponds to a distance from the 
Sun of 2754 pc for NGC 6819, 1259 pc for NGC 2099, 912 
pc for NGC 2168, and 1096 pc for NGC 2323. The surveyed 
areas are then 182, 81, 56 and 72 pc^, respectively. The cor- 
responding astrated mass per unit area is then 23.1, 32.1, 
15.2 and 6.7 Mqpc~^, equivalent to star formation rates of 
4.6x 10"^ 6.4x 10"^ 3.0x 10"^ and 1.3 x lO-'^Afoyr" V^ 
respectively. These values refer only to stars with masses in- 
cluded in the adopted sets of evolutionary tracks, namely 
from 0.6 to 100 Mq. To derive total astrated masses and 
star formation rates, one should extrapolate down to 0.1 
Mq, following the preferred IMF. In practice, extrapolating 
Salpeter's IMF down to 0.1 Mq, we must multiply by a fac- 
tor 2.05 the astrated masses and the rates. If, at the other 
extreme, we prefer to assume the positive slope -1-0.44 below 
0.6 Mq suggested by Gould, BahcaU & Flynn (1997), then 
the above values must be multiplied only by 1.27. 

There is a significant difference (more than a factor of 
four) in the star formation rates of our four clusters, the 
most active having been NGC 2099 and the least active 
NGC 2323, with no apparent relation with Galactic loca- 
tion. These rates are 2-3 orders of magnitude higher than 
the average star formation rate of field stars in the Solar 
neighbourhood, (2 — 10) x 10~^ MQyr~^pc~^ (e.g. Timmes, 
Woosley & Weaver 1995), and are comparable to those of 
the most active dwarf irregular and blue compact galaxies 
(e.g. Tosi 2003 and references therein). The survival factor of 
the stars formed in any cluster clearly depends both on age 
and on evaporation and NGC 6819 turns out to be the only 
one of the four clusters having suffered a reduction from the 
number of formed stars to that of the stars still alive and 
present in the studied area. 



7 CONCLUSIONS 

We re-derive key parameters for the four very rich open star 
clusters, NGC 6819, NGC 2099, NGC 2168, and NGC 2323. 
The parameters are measured by comparing high quality. 
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Table 2. Summary of Parameters and Results 



Parameter 



Explanation 



NGC 6819 



NGC 2099 



NGC 2168 



NGC 2323 



Position: 

c«J2000 

'5j2000 

^J2000 

6j2000 



Distances: 

{m—M)Y - Previous Study 
{m-M)v - FST (C = 0.03) 
(m-M)v - BBC94 
(m-M)v - FRANEC 



RA 

declination 
Galactic longitude 
Galactic latitude 



apparent distance modulus 
apparent distance modulus 
apparent distance modulus 
apparent distance modulus 



19'^41™17.7= 
+40° 11' 17" 
73.98° 
8.47° 



12.30 ± 0.12 
12.31 
12.57 
12.47 



06'^08™54.0= 
+24°20.0' 
177.65° 
3.09° 



11.55 ± 0.13 
11.52 
11.62 
11.42 



06'^08™54.0= 
+24°20.0' 
186.58° 
2.18° 



10.42 ± 0.13 
10.22 
10.42 
10.48 



071^ 02™ 48.0= 
-08°22'36" 
221.67° 
-1.24° 



10.68 ± 0.14 
10.68 
10.88 
10.68 



E{B-V) 
E{B-V) 
E{B-V) 
E{B-V) 



Previous Study 
FST (C = 0.03) 
BBC94 
FRANEC 



reddening 
reddening 
reddening 
reddening 



0.10 
0.10 
0.12 
0.15 



0.21 
0.36 
0.36 
0.36 



0.20 
0.20 
0.20 
0.22 



0.22 
0.22 
0.22 
0.22 



(m-M)o 
(m-M)o 
(m-iW)o 
(m-M)o 



Previous Study 
FST (C = 0.03) 
BBC94 
FRANEC 



d - Previous Study 
d - FST (C = 0.03) 
d - BBC94 
d - FRANEC 



true distance modulus 
true distance modulus 
true distance modulus 
true distance modulus 

distance from Sun 
distance from Sun 
distance from Sun 
distance from Sun 



11.99 ± 0.18 

12.0 

12.2 

12.0 



2500 ± 
2512 pc 
2754 pc 
2512 pc 



216 

199 V'- 



10.90 ± 0.16 
10.40 
10.50 
10.30 

1513 ± III pc 
1202 pc 
1259 pc 
1148 pc 



9.80 ± 0.16 
9.60 
9.80 
9.80 



912 ± 
832 
912 
912 



pc 



10.00 ± 0.17 
10.00 
10.20 
10.00 



1000 ± 
1000 
1096 
1000 



81 



pc 



Age: 

t - Previous Study 
t - FST (C = 0.03) 
t - BBC94 
t - FRANEC 



isochrone fit age 
synthetic CMD age 
synthetic CMD age 
synthetic CMD age 



2.5 Gyrs 
2.3 Gyrs 
2.0 Gyrs 
1.6-1.8 Gyrs 



520 Myrs 
520 Myrs 
430 Myrs 
400 Myrs 



180 Myrs 
200 Myrs 
180 Myrs 
110 Myrs 



130 Myrs 
130 Myrs 
120 Myrs 
110 Myrs 



Metallicity: 
Z 



heavy metal abundance^ 0.02 



<0.02 



0.012 



0.020 



Binary Fraction 
Binaries 



binary percentage 



20% 



20% 



30% 



30% 



Size: 

e 

D 
M 



angular diameter 
linear diameter'^ 
mass of cluster'' 



19' 

15.2 pc 
4000 Mq 



27.'8 
10.2 pc 
2500 Mq 



>32' 
>8.5 pc 
850 Mq 



30' 

9.6 pc 

480 Mq 



1. These values were not spectroscopically determined and reflect those used in the original isochrone fits. 
2. Computed using the distances from the BBC94 results. 
3. M is the mass of gas that formed stars between 0.6 and 100 Mq and is therefore only a lower limit to the total astrated mass. 



deep empirical colour-magnitude diagrams with MonteCarlo 
simulations. The combination of comparing the morphol- 
ogy and number of stars in various evolutionary phases and 
cluster luminosity functions allows us to provide tight con- 
straints on the reddening, distance, age and binary fraction 
for each cluster. In cases where the cluster metallicity is not 
certain, simulations with different abundances are also com- 
pared. In all cases the data are better reproduced when a 
fraction of unresolved binary systems between 20 and 30% 
is assumed. A summary of the results is given in Table 2. 

The synthetic CMDs and LFs are generally found to be 
in excellent agreement with the observational data. This cir- 
cumstance, combined with the fact that different sets of stel- 



lar evolution tracks provide different values for the cluster 
parameters, confirms how important it is to use more than 
one set of models to estimate the theoretical uncertainties. It 
also shows that a homogeneous approach is crucial to derive 
reliable overall cluster properties, such as age-metallicity re- 
lations. In fact, cluster dating based on different stellar mod- 
els may lead not only to different absolute ages, but also to 
different age ranking, with significant drawbacks on the in- 
terpretation of the cluster properties in terms of Galactic 
evolution. 

Some discrepancies, such as the thickness of the main 
sequences vis a vis the size of the photometric errors, and 
the possibility that cluster stars may be contaminating our 
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blank field (particularly in the case of NGC 2168) are iden- 
tified and discussed. 

For each cluster, we also measure the astration mass 
(i.e., the total mass that went in all of the stars formed 
in the cluster) according to a single slope Salpeter IMF. 
From this, we calculate the star formation rate between 0.1- 
100 Mq to be 9.4 x 10"®, 1.3 x 10-^ 6.2 x 10"^ and 2.7 x 
lO-'^Moyr-V"^ for NGC 6819, NGC 2099, NGC 2168, 
and NGC 2323 respectively. These rates would be a factor 
1.6 lower if the IMF below 0.6 Mq has a slope of +0.44, 
as inferred by Gould et al. (1997) from HST data, rather 
than Salpeter. The true value may in fact lie somewhere 
in between these two IMFs, as recently discussed e.g., by 
Chabrier (2003). In principle our data are deep enough to 
allow for a direct derivation of the cluster IMF; however, the 
four clusters are heavily contaminated by fore/background 
stars and the lack of appropriate decontamination fields has 
prevented us from a safe analysis of the star counts at the 
fainter magnitudes. The need of appropriate photometry in 
nearby fields is emphasised, also for the purpose of adequate 
studies of the clusters mass segregation/evaporation. 
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